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We investigate the zero-temperature phase diagram of the nearest-neighbor kagome antiferro- 
magnet in the presence of Dzyaloshinksii-Moriya interaction. We develop a theory for the transition 
between Z2 spin liquids with bosonic spinons and a phase with antiferromagnetic long-range order. 
Connections to recent numerical studies and experiments are discussed. 



I. INTRODUCTION 

The nearest neighbor spin S = | antiferromagnet on 
the kagome lattice has been the focus of extensive the- 
oretical and experimental study because it is a prime 
candidate for realizing a ground state without antiferro- 
magnetic order. 

On the experimental side, much attention has fo- 
cused on the S = 1/2 compound herbertsmithite 
ZnCu 3 (OH) 6 Cl 2 . It has J « 170 K and no observed 
ordering or structural distortionP^ However, there is an 
appreciable amount of substitutional disorder between 
the Zn and Cu sites (believed to be of the order 5-10 %) 
which affects the low T behaviorPE*! More importantly, 
there is an upturn in the susceptibility at T = 75 K which 
has been ascribed to the DM interaction d 5 * 12 * ^. 

On the theoretical side, the most recent evidence^HH] 
on the the nearest-neighbor antiferromagnet points con- 
sistently to a ground state with a spin gap of 0.05 J and 
valence bond solid (VBS) order. The pattern of the VBS 
order is quite complex with a large unit cell, but was 
anticipated by Marston and Zeng22l by an application of 
the VBS selection mechanism described in the 1/N ex- 
pansion of the SU(iV) antiferromagnet 

The influence of the D M interactions has also been 
studied theoretically22E3H2D Starting with an "algebraic 
spin liquid" ground state, Hermele et al^ argued that 
the DM coupling, D, was a relevant perturbation, im- 
plying that an infinitesimal D would induce long-range 
magnetic order. In a recent exact diagonalization study, 
Cepas et al^ reach a different conclusion: they claim 
that there is a non-zero critical DM coupling D c be- 
yond which magnetic order is induced. They estimate 
D c /J « 0.1, quite close to the value measured^ for 
ZnCu 3 (OH) 6 Cl 2 which has D/J w 0.08. This proximity 
led Cepas et al. to suggest that the quantum criticality 
of the DM-induced transition to magnetic order controls 
the observable properties of this kagome antiferromagent. 

The purpose of this paper is to propose a theory for the 
quantum critical point discovered by Cepas et aZP3. We 
will compute various observables of this theory, allowing 
a potential comparison with numerics and experiments. 

Given the evidence for VBS order in the model without 
DM interactions^^, it would appear we need a theory 
for the transition from the VBS state to the magnetically 



ordered state. However, the VBS ordering is weak, and 
can reasonably be viewed as a perturbation on some un- 
derlying spin -liqui d ground state. Schwandt, Mambrini, 
and Poilblan d 20 * 21 ! have recently presented evidence that 
the kagome antiferromagnet is proximate to a Z2 spin 
liquid state, and that vison condensation in this state 
leads to weak VBS ordering. Their dimer representation 
leads naturally to Z 2 spin liquid states in the same class 
as that orig inally describe d 27 * 28 ! by the Schwinger boson 
methocP^ESI. We will therefore neglect the complexities 
associated with the VBS ordering and work with the par- 
ent Zi spin liquid state. This is equivalent to ignoring 
the physics of the vison sector, and assumes that the 
magnetic ordering transition can be described in a the- 
ory of the spinons alone^. The main result of this paper 
will be a theory of the quantum phase transition from 
the Schwinger boson Zi spin liquid to the magnetically 
ordered state as induced by the DM interactions. 

We will begin in Section [II] with a description of the 
mean-field theory of the Z 2 spin liquid and its transi- 
tion to the magnetically ordered state in the presence of 
DM interactions. This will be ca rried using the Sp(Y) 
Schwinger boson formulation^*^!, for which the mean- 
field theory becomes exact in the large N limit. We 
will turn to fluctuation corrections and the nature of the 
quantum critical point in Section |III| Here we will show 
that the critical theory is the familiar three dimensional 
XY model. However, its connection to experimental ob- 
servables is subtle: in particular, the XY field itself is not 
directly observable. 

While this paper was in preparation, a description of 
the Schwinger boson mean field theory in the presence of 
DM interactions also appeared in Ref. they consider 
mean-field solutions with larger unit cells than we do, 
but did not analyze the critical field theory. Where they 
overlap, our results are in agreement with theirs. We 
also note the recent experimental observations of Helton 
et aZP% who present evidence for quantum criticality in 
ZnCu 3 (OH) 6 Cl 2 . 



II. MEAN FIELD THEORY 

The model we consider is a standard Heisenberg Hamil- 
tonian supplemented by an additional DM interaction. It 
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FIG. 1: Staggered DM interaction from triangle to triangle 
in z-direction as in Ref. 1341 The arrowheads indicate D to 
come out of the plane, whereas the tails denote D to go into 
the plane. Note that the triangles are summed clock- and 
anticlockwise, respectively, indicated by the arrows on the 
bonds. 



assumes the form 



(1) 



Si in this notation denotes the spin operator at site i, 
Jij is assumed to be uniform and of the nearest neigh- 
bor type, and D,j = Dije z is taken along the z-axis and 
staggered from triangle to triangle^, see Fig. [TJ This 
additional term explicitly breaks the spin-rotation sym- 
metry by favoring configurations lying in the x-y plane. 
Furthermore this term increases the tendency o f clas sical 
spin ordering. It has been shown in earlier work a 29 l 3 ° lthat 
Schwinger bosons are ideally suited to describe phase 
transitions between paramagnetic and magnetically or- 
dered phases in spin models. Following Ref. [37|we intro- 
duce a Sp(iV) generalization of the spin operators, which 
formally allows to consider a controlled large N limit par- 
ticularly suited for the study of frustrated spin systems 
such as kagome or traingular antiferromagnets. In the 
Sp(l) case, which is isomorphic to SU(2), one can repre- 
sent the spin variables as 



S. = b*„-^-b,„, 



with f being the Pauli matrices and with 



and 



J2 b 



25 = rib 



(2) 



(3) 



In the case of the large- N generalization the Schwinger 
bosons acquire another index counting the copies of the 
system (we drop this index in the following discussions, 
but display N whenever it is essential). The large- N 
generalization formally justifies the mean- field, with the 
saddle point becoming exact in the limit N — > oo. In 
order to properly reformulate the problem at hand we 



introduce two decoupling parameters 



Qij ^ ^ ^aa'bio-bju' 

era' 

Pij = y ] Tgg'bicrbja-i 



(4) 



where e CT<T < is the antisymmetric tensor and t x is just the 
standard Pauli matrix. We see from the above expres- 
sions that Qij = —Qji whereas Pij — Pji. This implies 
that the bond variables do not have a direction. 

The constraint Eq. ([3| is implemented via a La- 
grangian multiplier in a standard way. The Hamiltonian 
of the system formulated in the fields defined in Eq. Q 
consequently reads 



H 

N 



(5) 



where K = n^/N. We furthermore introduce N u as the 
number of unit cells in the systems and N s as the number 
of sites within the unit cell. We can write the mean-field 
Hamiltonian per flavor and unit cell as 



H 



AIF 



^^**(k)H( %i , Pii ,A,k)*(k) 



J 
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N S \{1 + k) 



iD 



(ij) 1 



(6) 



where Yluj)' denotes the sum over bonds belonging to 
the unit cell, 

**(k) = (^ T (k),...,^ st (k),6 u (-k),...,6^(-k)) , (7) 
and the matrix 

/ AI C*(fc) 



^ C(Jfc) AI 



(8) 



with the matrices I (identity) and C(A;) being N s x N s 
matrices; the explicit form of these matrices is given in 
Appendix \K\ As mentioned before, one of the major as- 
sets of the Schwinger boson approach is that it can de- 
scribe magnetically disordered gapped spin liquid phases 
as well as magnetically ordered states. On a formal level 
in the large N approach this is achieved by introducing 
the following notation for the Schwinger bosons 



(9) 



(VNxi^Y 1 ) where m = 2,...,N. 



The first component is thus a classical field. If x ^ 
it signals condensation which causes long range order to 
appear. Following Ref. 1271291 we can integrate out the 
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O u 
□ v 
O w 



This states is therefore called the k = state. 
(ii) qi = —q2 real, p\ = —p^ pure imaginary: upon in- 
creasing k, the Schwinger bosons condense at wavevector 
k = ±(27r/3, 0) into a state which is called the -\/3 x \/3 
antiferromagnet, characterized by an enlarged unit cell. 
Solutions with larger unit cells can be present with addi- 
tional frustrating interactions^, but we will not consider 
them here. 




FIG. 2: We take a unit cell of three sites, labeled u, v, w. 
The arrows indicate the values of the oriented variables gy. 
The links with single arrows have qy — qi, while those with 
double arrows have gy = q2- The Py are unoriented, and 
take values p\ and p2 on these links respectively. 



Schwinger bosons and the zero-temperature mean field 
energy assumes the following form 



k,H=l,..,N. 



J 



+ o 0®i| 2 - (Qij^a'XiaXja') + h.C.)] 



iD 



^(PijQii ~ QijPij) 



(«)' 



iD 



^ ) iPij^cru'XjcrXjui + qijT cra ,X i(J Xj (J i) 



(«)' 



iD 



(10) 



where denotes the sum over all sites within one unit 
cell. In the following we solve the self-consistency equa- 
tions according to 



K = ( b ia b ia)MF 
Qij = (Qij)MF , 
Pij = (Pij)MF ■ 



(11) 



with the Hamiltonian defined in Eq. (J6|. 

Our s olutio n of the mean-field equations follows previ- 
ous workp212Sl which classified physically different Z2 spin 
liquid solutions without the DM interactions. We found 
that these solutions have a natural generalization in the 
presence of DM terms, with values of the pij which re- 
flect the symmetries of the q^ . Two stable solutions were 
found in previous work, with only two possibly distinct 
values of q^ as illustrated in Fig. [2] Including the DM 
interactions, these solutions extended to 
(i) qi — <72 real, p\ — P2 pure imaginary: upon increas- 
ing k, the Schwinger bosons condense at k = 0, with the 
spins at angles of 120° to each other within the unit cell. 



A. Phase diagram 



Our phase diagram is shown in Fig. |II A| as a function 
of k = rib/N (which corresponds to the spin size) and the 
parameter D / J . Our phase diagram is similar to that ob- 

D/J 




FIG. 3: Phase diagram in the iV — > 00 limit of the Sp(iV) 
theory. The x-axis shows k = 2S/N and the y-axis D/J. The 
phases have long-range magnetic order type (LRO), or gapped 
Z2 spin liquids. The thick line is a first-order transition, while 
the thin lines are second-order transitions: the transition to 
the k = LRO state is described in Section ULTl The limit of 
D I J — > reduces to the results presented in Ref. [27| 

tained recently by Messio et alW^. They also considered 
solutions with a larger unit cell which were stable over 
some portion of the phase diagram. 

We start with a discussion of the classical limit. While 
for D = the long-range ordered state of the -\/3~ x 
type is generically preferred^, infinitesimal D favors the 
so-called k = state. This is reproduced by our mean- 
field equations in the large spin limit k — > 00. For finite 
values of k there is a small slab in which long-range order 
of the \/3 x \/2> type is favored over the k = 0. These 
states are separated by a first order transition driven by 
the ratio D/J. 

A similar behavior appears for the corresponding spin 
liquid states at small k. The q\ = —q2 is favored at small 
D/J, and then undergoes a first order transition to the 
Qi = 92 state at large D/J. 

In exact diagonalization studies of the spin S = 5 
kagome antiferromagnet with DM interactions a second 
order phase transition between a phase with short ranged 
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k = correlation d 18 ! 35 ! (obtained in the pure Heisenberg 
case) and phase with k = long range order was founcpH 
Such a transition is also present in our mean-field theory, 
which therefore can be used for a study of critical prop- 
erties in Section UTIl 



III. QUANTUM CRITIC ALITY 

We will consider only the transition out of the q\ = q 2 
spin liquid, because that is what is seen in the numer- 
ical studies^!. The corresponding transition out of the 
qi = —q 2 spin liquid can be treated in a similar man- 
ner. Throughout this section we will consider the physi- 
cal SU(2) antiferromagnet directly, and not take the large 
N limit. The method followed below has been reviewed 
in a more general context in Ref. 1301 

Since we are at k = 0, we can write the effective action 
for the bosons by making the small momentum expansion 
of the matrix in Eq. (A2). We take 3 sites, u, v, w, in 
each unit cell (see Fig. [2J), and then take the continuum 
limit of the saddle-point Lagrangian. We write the boson 
operators on these sites as b ua — U a , b va — V a etc., and 
set qi = q2 = q, and pi — p 2 = ip with q and p real. 
Then, we write the final Lagrangian in the form 



£ = £ 



H 



C 



DM 



(12) 



representing the contributions of the Heisenberg ex- 
change and the DM coupling, respectively. 

From Eq. ( A2 ) we obtain the Lagrangian in the ab- 



sence of a DM term (which describes the k = solution 
of Ref. OH]): 



dU n 



v: 



dr 



dr 

-x(\u a \ 2 + \v a 

- Jqe aa > {U a V a > 



W, 



Or 
f |W CT | 2 ) 
V a W a . 4 



Jq 



W a U a .) +c.c. 
{dxUvdxVv- + d 2 V a d 2 W a , + d 3 W a d 3 U a ,) 

(13) 



+ c.c. 



where di is the gradient along the direction e.; in Eq. ( |Al| . 

We now perform a unitary transformation to new vari- 
ables X a , Y a , Z n . These are chosen to diagonalize only 



the non-gradient terms in Ch- 




(14) 



where £ = e 27 ™/ 3 . The tensor structure above makes it 
clear that this transformation is rotationally invariant, 



and that X„, Y„, Z a transform as spinors under SU(2) 
spin rotations. Inserting Eq. (14) into Ch we find 



£ 



H 



X. 



dX n 



y:^ + z:^ + (x + vsjq) i z a i 



a dr a dr 
I- (A - V3Jq)\Y a \ 2 

+ {\dxZa 



■ x\x a 



2 + \dyZ„\ 2 ) 



(15) 



The ellipses indicate omitted terms involving spatial gra- 
dients in the X a and Y a which we will not keep track of. 
This is because the fields Y a and X a are massive rela- 
tive to Z a (for q < which is the case in our mean-field 
solution), and so can be integrated out. This yields the 
effective Lagrangian 



C z 



1 



H 



-\d T Z a \ 



JqV3 



(A - Jq) 

+ (X + V3Jq)\Z a \ 2 



(\d X Z„\ 2 + \dyZ a \ 



(16) 



Note that the omitted spatial gradient terms in X a , Y a 
do contribute a correction to the spatial gradient term 
in Eq. (16), and we have not accounted for this. This 



Lagrangian shows that the mean-field theory has a tran- 
sition to magnetic order at A = |\/3Jq|, which agrees 
with earlier resultJ^H. 

The effective Lagrangian £ H is almost the complete 
solution for the critical theory in the system without the 
DM interactions. However, we also need higher order 
terms in Eq. (16), which will arise from including the 
fluctuations of the gapped fields Q and A. Rather than 
computing these from the microscopic Lagrangian, it is 
more efficient to deduce their structure from symmetry 
considerations. The representation in Eq. (14), and the 
connection of the U, V, W to the lattice degrees of free- 
dom, allow us to deduce the following symmetry trans- 
formations of the X, Y, Z: 

• Under a global spin rotation by the SU(2) matrix 
gaa', we have Z a ga-a-iZa-/, and similarly for Y, 
and Z. When DM interactions are included, the 
global symmetry is reduced to U(l) rotations about 
the z axis, under which 



i0 r 



Y t -> e lU Y t 
X t e w X t 



Y i -> 

X i -> e-*%. 



(17) 



• Under a 120° lattice rotation, we have U a 
V a W a -> U a . From (T2J), we see that this 

symmetry is realized by 



Z„ 



CZa ! Ya CXa , X a 



x„ 



(18) 



Note that this is distinct from the SU(2) rotation 
because det(C) ^ 1. 
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Under time-reversal, we have U a — » t aa iU*,, and 



similarly for V a , W a . This is realized in Eq. (14) 

by 



Z a -t iZ a , Y„ -> —iY a , X a -> e a<J >X* 



(19) 



In particular, note that does not map to Z^ 
under time-reversal. 



It is easy to verify that Eq. (151 is invariant under all the 



symmetry operations above. These symmetry operators 
make it clear that the only allowed quartic term for the 
Heisenberg Hamiltonian is (52 a \Z a \ 2 ') ■ this implies that 
the Z2 spin liquid to antiferromagnetic order transition of 
this model is in the universality class of the 0(4) modeP^. 
Let us now include the DM interactions. From 



Eq. (A2), we see that 



C DM = i*yTZ al (U a V a , + V a W a , + W^+c.c. (20) 
We have dropped a term proportional to Dp which has 



the same structure as the terms in Eq. (151, and ignored 



spatial gradients. In terms of the fields in Eq. (14 1, this 
takes the simple form 



C-DM = ~Y {ir^ a ,X a X a i + C.C. 



(21) 



the energy densities and does not have any important ef- 
fects. More interesting is the w term, which shows that 
the global symmetry is not 0(2)<g>0(2) but 0(2)<g)Z 12 . In 
the magnetically ordered phase with (Zf) ^ (say), this 
term will induce a small ordering field ~ Z® in the XY 
model for Z±. However, the action for Zi has a 'mass' 



term with a positive co-efficient, and this sixth order 
term will not immediately induce ordering in Z^\ i.e. a 
magnetic phase with (Zf) 7^ and (Z±) = has a fi- 
nite range of stability. Thus close to the transition we 
can neglect the Z± field entirely, and transition is in the 
universality class of the three-dimensional XY model. 

The choice above of Z-f over Z± gives the incorrect 
appearance that we are breaking the spin reflection sym- 
metry S z —J- —S z of "H, suggesting the appearance of a 
net z ferromagnetic moment. However, notice that the 
theory of is relativistic, and so contains both spinons 
and anti-spinons which carry S z = +1/2 and S z = —1/2 
respectively. The spinon of Z^ also carries S z = —1/2, 
and this is degenerate with the anti-spinon of Z^ in the 
0(4) invariant theory in Eq. (16). It is this latter de- 
generacy which is lifted by the DM interactions, which 
induce a vector spin chirality along the z direction^ (as 
we will see below). We will also see there is no net fer- 
romagnetic moment, because time-reversal symmetry is 
preserved. 



Observables 



and it can be verified that these terms are invariant under 
Eqs. [ |17|18|19p . As before, we now integrate out X a and 
Y„ from Ch + C-dm in Eqs. (15) and (21). We obtain 
a Lagrangian with the same structure as Eq. (16), but 



all couplings become dependent upon <r; in other words, 
we have 2 separate XY models for Z^ and Z±. Perform- 
ing a careful analysis of allowed higher order terms as 
restricted by the symmetry constraints discussed above, 
and after appropriate rescalings of the spatial, tempo- 
ral, and field scales, we obtain the field theory with the 
Lagrangian 



C z = 



\9rZ t \ 2 
\drZ.l 2 



v\z,\ 2 \za 2 



VZ t \ 2 + s t \Z t \ 2 + u t |Z t | 4 
VZ^ + s^Ztf + u^f 

w^Z^f + iZ^Zlf). (22) 



Note that s-j- 7^ sj, in general (and similarly for u^^ etc.), 
and equality is not required by the time-reversal symme- 
try in Eq (19). Time-reversal symmetry does prohibit a 
term ~ (ZfZ±) which is allowed by the other symme- 
tries. Thus we expect only one of Z+ or Z^ to condense at 
the quantum critical point: as we will see from the anal- 
ysis of observables in Section III A this transition does 



indeed correspond to the development of spiral magnetic 
order in the x-y plane. The choice between Z^ and Z^ is 
controlled by the sign of D. 



Eq. ( 22 ) also contains terms which couple the two XY 



To determine the operators corresponding to the ferro- 
magnetic moment, let us couple a uniform external field 
h to the lattice Hamiltonian. This adds to the continuum 
Lagrangian the term 



-h ■ Taa , (I/* u a . + v; v a , + W*W a , ) (23) 



Inserting the parameterization in Eq. ( 14 ) this becomes 

c h = -h • r aa , (x;x a , + y;z a , + z*Y a ,) (24) 



We now need to integrate out X a and Y a in the La- 
grangian Ch+C-dm+Ch, defined by the sum of Eqs. (15), 



(21 ), and (24), and collect the terms linear in h. Without 



the DM coupling, we obtain 



dr 



(25) 



models to each other. The lowest allowed term, v, couples 



Comparing with Cfj in Eq. (Il6|) we see that this is just 
the coupling to the conserved~~SU(2) charges of the 0(4) 
model: this is the usual term which determines the mag- 
netic susceptibility of the Heisenberg antiferromagnetP3. 
Upon including the effects of Com we find that the 
essential structure of Eq. (25) does not change: the 
r matrices get multiplied by some tr-dependent factors 
T aa' — > fo-Taa' fa' which do not modify the scaling con- 
siderations. No term with a new structure is generated 
by the DM coupling. It can now be seen that these ex- 
pressions have vanishing expectation values under Cz in 



G 



Eq. ( 22 ) , and so there is no net ferromagnetic moment in 



the absence of an external field. 

We now turn to the antiferromagnetic order parameter; 
for a coplanar antiferromagnet, this is described by 



S t oc Ni cos(Q • Ti) + N 2 sin(Q • r 4 ), 



(26) 



where Ni j2 are 2 orthogonal vectors representing the spi- 
ral order, and Q is wavevector at which the spin structure 
factor is peaked. For our model, we can see that 



(27) 



Using Eq. ( 14 ) , and keeping only the lowest order term, 
we therefore obtairP-2 



the physically measurable magnetic order was the square 
of the XY field. 

The behavior of the uniform magnetic susceptibility 
follows from the scaling dimension of the operators in 
Eq. (25). For h along the z direction, the magnetization 
is the just the conserved U(l) charge of the XY model: 
scEl it has scaling dimension 2, and the susceptibility ~ 
T. For h along the x or y directions, we have to integrate 
out Zj,, and then the lowest dimension operator coupling 
to the square of the field is |Z^| 2 . This means that the 
susceptibility only has a weak singularity at the quantum 
critical point given by that in (|Zj-| 2 ): at the quantum 
critical point, there is a non-analytic term ~ j>3— l/i^ 
above an analytic background. 




n, = ! -(z 2 + z T2 )/2 



(28) 



in a notation that makes the rotational invariance evi- 
dent, this relationship is 



Ni + iN 2 = -e a/ 3Tp a Z a Z a 



(29) 



Note that a phase with (Zf) ^ and (Zj,) = has spiral 
order in the x-y plane. 

To complete the list of operators which are quadratic 
in the Z CT! we consider the vector spin chirality^. This 
is defined here by 



S u x S v 4- x S w S w x S i , 



(30) 



Using Eq. ( 14 ) we find that the leading operator mapping 



to vector spin chirality is (dropping an overall factor of 
l^tl 2 + l^l 2 ) 



Z* a T Ga i Z a i . 



(31) 



Note that in the presence of the DM term, the couplings 
in the effective theory (22 1 imply that the z component 



of the vector spin chirality is always non-zero. 



IV. CONCLUSION 

We have presented a theory for the quantum critical 
point between a Z 2 spin liquid and an ordered antiferro- 
magnet for the kagome antiferromagnet in the presence of 
DM interactions. The critical theory is just the three di- 
mensional XY model. However, the XY order parameter 
carries a Z 2 gauge charge, and so it is not directly observ- 
able. In particular, the antiferromagnetic order parame- 
ter is the square of the XY order parameter. Specifically, 
the theory is given by Lz in Eq. p2| , and its observables 
are described in Section fill Al 

It is interesting to compare our results with recent ob- 
servations of quantum critical scaling in ZnCu3(OH)6Cl 2 
by Helton et al. 33 . Their neutron scattering measure- 
ments show an antiferromagnetic susceptibility which 
scales as T -0 - 66 . This is actually in reasonable agreement 
with our theory, which has a susceptibility ~ j"-o.526 
However, they also observe a similar divergence in mea- 
surements of the uniform magnetization, while our theory 
only predicts a very weak singularity. We suspe ct that 
this difference is due to the present of impuritieJ ^ 9 * 15 * 4 ^ , 
which can mix the uniform and staggered components. A 
complete study of impurities near the quantum critical 
point described above is clearly called for. 



B. Critical properties 

Let us assume the transition to magnetic order pro- 
ceeds via the condensation of Zf. The transition is in 
the XY universality class, and the dimension of the an- 
tiferromagnetic order parameter is 

dim[N x ] = dim[N 2 ] = dim[Z|] = . (32) 

The value of the exponent rj can be rea d off from re- 
sults for the three-dimensional XY modeP^H, and we 
obtain rj as 1.474. The antiferromagnetic susceptibility 
will diverge at the critical point as T~ < - 2 ~' , ' ) . We note the 
recent work of Ref. HU] in a different context, which also 
considered a model with an XY critical point at which 
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Appendix A: Microscopic form of the Hamiltonian 

We here give explicit expressions for the Hamiltonian 
introduced in Sec. [Hj We introduce the following set of 
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unit vectors 

e 3 = o(-l,0) (Al) 



Appendix B: Dispersion of the lowest excited spinon 

state 



For qi — q2 and D — 0, the ground state is doubly 
degenerate. The degeneracy splits as D moves away from 
0. In all cases, the minimum excitation energy occurs at 
k = 0. 

Fig. [4] plots the momentum dependence of the energy of 
the lowest excited spinon at D/J = 0.3, k = 0.2. There 
is a finite energy gap at k = making it a gapped spin 
liquid. On the other hand, Fig. [5] is the dispersion plot 
for a case with long-range magnetic order. The energy 
gap at k = is closed and condensation occurs at this 
wave vector. 



For qi = —q-2, there is a unique ground state even for 
D = 0. Energy minima is at k = ±(2tt/3,0). Fig. [6] 
shows the dispersion of the lowest lying state of the spin 
liquid with D/J — 0.03, k = 0.55. A case with long range 
ordering is shown in Fig. [7J 



which allows to express ki = k • e^. In the following 
we concentrate on the states with 3 sites per unit cell. 
In that case the vector \E' introduced in Eq. (J7]) has 6 
components and the matrix C consequently is a 3 x 3 
matrix with the following entries: 



(A2) 
I 




4 



FIG. 4: Momentum dependence of the energy us(k) of the 
lowest excited spinon state of the kagome-lattice quantum 
antiferromagnet for the q\ = q2 state at D/J = 0.3, k = 0.2. 
The minimum excitation energy is at k = and has a finite 
energy gap. 
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FIG. 5: Dispersion of the lowest excited spinon state for the 
qi = q2 state at D/J — 0.3, k = 0.4. The energy gap closes 
at k = and condensation occurs. 




FIG. 6: disorder: q\ = -q 2 , D/J = 0.05, k = 0.3 



Appendix C: Condensation 



tion is energetically favored while the opposite is true for 
D < 0. The two condensations are degenerate for D = 0. 

For the q± = —q-2 state, condensation occurs at k = 
k = (2ir/3, 0) or k — —k. The two states have identi- 
cal energies and the condensation is spontaneously cho- 
sen. Similar analysis to Eqns. (13) and (14) gives the 




FIG. 7: order: q ± = -q 2 , D/J = 0.03, n = 0.55 



eigenvectors corresponding to the lowest lying state. For 
k = k this is (i, —i, i, 1, —1, 1) while for k = — k the cor- 
responding eigenvector is (— — i, 1, — 1, 1). Therefore 
condensations can be parametrized as 
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(C2) 



For the qi = qi state Eq. (14) leads to the following 
parametrization of the condensation of Zj. field at k = 0. 



1 

c 



(CI) 



where £ is the size of the condensate. Condensation of 
can be written similarly. For D > 0, field condensa- 



for k — k, and 



for k 
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